Magnetosonic solitons in a Fermionic quantum plasma 
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Starting from the governing equations for a quantum magnetoplasma including the 
quantum Bohm potential and electron spin-1/2 effects, we show that the system of 
quantum magnetohydrodynamic (QMHD) equations admit rarefactive solitons due 
to the balance between nonlinearities and quantum diffraction/tunneling effects. It is 
found that the electron spin-1/2 effect introduces a pressure-like term with negative 
sign in the QMHD equations, which modifies the shape of the solitary magnetosonic 
waves and makes them wider and shallower. Numerical simulations of the time- 
dependent system shows the development of rarefactive QMHD solitary waves that 
are modified by the spin effects. 

PACS numbers: 52.27.-h, 52.27.Gr, 67.57.Lm 



I. INTRODUCTION 



lere is currently a great deal of interest in collective quantum effects in plasmas 
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progress and techniques 
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these studies are motivated by recent experimental 
and also by possible astrophysical applications 
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Il9l . |20| In particular, magnetohydrodynamic (MHD) plasmas are of interest in 



such astrophysical applications. However, in strong magnetic fields, single electron effects 
that depend on the electron spin properties, such as Landau quantization, will be important. 
It is thus not surprising that collective s pin effects can influence the wave propagation in 
a strongly magnetized quantum plasma [l2j, [2l|, |22j, |23j. Moreover, the recent progress in 
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25| may offer an interesting 



producing ultra-cold plasmas in terms of Rydberg states 
experimental environment for quantum plasma dynamics. In such cold plasmas, the thermal 
energy of the particles can be very small compared to the Zeeman energy of the particles in 
an magnetic fields. Thus, collective spin properties of quantum plasmas may be possible to 
detect in a near future. 
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In this Brief Report, we will show that the balance between the nonlinear plasma and 
quantum effects gives rise to magnetosonic solitons. Using the governing equations for 
QMHD plasmas with tunneling and spin effects included, we derive a Sagdeev potential for 
the one-dimensional system. We show that in a magnetized quantum plasma, the electron 
spin-1/2 effect can strongly modify the amplitude and width of rarefactive solitons. 



II. GOVERNING EQUATIONS 

We begin by presenting the general governing equations for a quantum magnetoplasma 
in which the electron— 1/2 spin effect are included. We define the total mass density p 
[in, ii, + m.j7ij), the center-of-mass fluid flow velocity V = (m e n e v e + mifiiVi) / p, and the 
current density j = —en e v e + en^. Here m e (mj) is the electron (ion) mass, n e (n$) is 
the electron (ion) number density, v e (vi) is the electron (ion) fluid velocity, and e is the 
magnitude of the electron charge. From the general set of spin-fluid equations [12| the 
corresponding QMHD equations can be derived [23]. From these, we immediately obtain 
the continuity equation 

f t +v.( P v) = o. (1) 

Assuming the quasi-neutrality, i.e. n e ~ rii, the momentum conservation equation reads 

P ( |r + V ■ v) V = j x B - VP + F Q , (2) 



dt 

where P is the scalar pressure in the center-of-mass frame, the current is given by j = 
Pq 1 ^ x (B — p: M), M = (u^p/mi)taiih.(p, B B/k B T e )B is the plasma magnetization due 



to the electron spin, and [121. 123] 



5!^ V (jf 2 Vp) + — tanh ( ^1) 



is the quantum force due to collective tunneling and spin alignment. Here fi B — eh/2m e is 
the magnitude of the Bohr magneton, h is Planck constant divided by 2ir, and c is the speed 
of light in vacuum. The generalized Faraday law takes the form 
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where r\ is the plasma resistivity. 



III. SPIN SOLITONS 

Next, we assume that the magnetic field is along the z direction such that B = B(x,t)z, 
while we have the velocity V = V(x,t)x and the density p(x,t). With this, the governing 
equations reduce to 
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dV y dV_ 

dt dx 
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dx 2 



(7) 



Here Ac = c/u c = h/2m e c is the Compton wavelength, uj c is the Compton frequency, 
C s = [kB(T e + Tj)/mj] 1//2 is the sound speed, A = i]/p is the magnetic diffusivity, the last 
term in Eq. ([6]) is the spin force divided by rrii, and we have neglected the inertial term in 
the Faraday law (J7J). 

If the resistivity is weak, we may neglect the last term in the Faraday law ([7]), and obtain 
the frozen-in-field condition p = p b, where b = B/B , with the background values denoted 
by the zero index. Then, Eqs. and ([7]) form a closed system, taking the form 
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{In [cosh [eb)\ + efetanh (eb)} 



(8) 
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and 

where we have introduced the Alfven speed Ca = (-Bo// i oPo) 1//2 and the temperature nor- 
malized Zeeman energy e = fiBB Q /kBT e . 

We now normalize our variables as i — u ci t, x = (uid/C^x = (u pi /c)x (where ui pi = 
(noie 2 ' i 'e^vrti) 1 ! 2 is the ion plasma frequency), v = V/Ca, and c s = C s /Ca- We then obtain 



dv d fv 2 \ db 2 d ^ 
dt dx \ 2 / dx s dx 



- 2 V^h^T (4r7Tr I +«I^-{ln[oosh(e6)]+e6tanh(e6)} ) (10) 
\cu ce \uJcOx \y/b dx I dx 



with v% = k B T e /rriiC A = (l/e)(n B B /miC\), and 

§ + >) = o, ,n, 

where we, for simplicity, drop the bars on the normalized coordinates. 

Next, we assume that v and b are functions of £ = x — t^t, where Vq is a constant speed 
(normalized by Ca)- Then Eq. (fTTj) can be integrated as t> = t>o(l — 1/6), where we used 
the boundary conditions 6=1 and t> = at |£| = oo, and Eq. fflOl) can be integrated twice 
to obtain 

(f)W)=», (12) 
where Z = y/b and the Sagdeev potential [26| for our purposes reads 
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(13) 



In deriving 0131) we have used the condition ^(1) = 0. In Figs. 1 and 2, we have plotted 
the Sagdeev potential as well as the profiles of the corresponding solitary waves for different 
sets of parameters. The solitary waves have only sub-Alfvenic speeds and are characterized 
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by a localized depletion of the magnetic field and density. In Fig. 1, we see that the solitary 
waves increase their amplitudes for smaller speeds. In the limit of zero speed, we have 
rarefactive solitons with a zero density at its center. The influence of the electron spin-1/2 
effect on the solitary waves is displayed in Fig. 2, where we see that larger values of e lead 
to wider solitary waves with shallower density and magnetic field depletions. In order to 
study the influence of the spin pressure on the nonlinear dynamics of our system, we have 
solved the time-dependent system of equations ([TO]) and (fTTl) for different values of the spin 
pressure parameter e. As an initial condition at t = 0, we took a magnetic field with a local 
depletion in the form Gaussian pulse 6 = 1 — 0.5 exp(— x 2 /100), while the velocity v was 
set to zero. For e = 5, we see in the left-hand column of panels in Fig. 3 that the initial 
pulse develops into two counter-propagating pairs of rarefactive solitary waves, where the 
smaller pulse in the pair propagates with a somewhat larger speed, ~ 0.75 Ca, than the 
larger one that propagates with a speed of ~ 0.65 Ca- For a larger value e = 10, displayed 
in the right-hand panels of Fig. 3, the pulse develops into two counter-propagating pulses 
that propagate with somewhat lower speed, ~ 0.4 Ca, and they are wider and of smaller 
amplitude than the large-amplitude pulses for e = 5. All pulses are rarefactive and are 
propagating with sub-Alfvenic speed, in agreement with our analysis in Figs. 1 and 2. 



IV. SUMMARY AND DISCUSSION 



In the numerical examples of the previous section, the normalized Zeeman energy e played 
a crucial role. In particular, the spin contribution to the soliton dynamics is enhanced when 
the Zeeman energy is of the order of or greater than one (we note however that other pa- 
rameters play a role in forming the necessary shape of the Sagdeev potential). Thus, it is 
natural to investigate what type of parameter values correspond to e > 1 . For astrophysical 
plasmas, such as in pulsar magnetospheres, we can have B Q < 10 10 T |20(], implying that the 
that e > 1 for T e < 10 9 K, i.e., not a very severe constraint. However, in such environments, 
the plasma often has relativistic temperatures and flows, and a relativistic formalism should 



be used. In the case of Rydberg plasmas 
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251 ]. where the temperature can go as low 



as millikelvins, we see that the Zeeman energy is greater than one for external magnetic 
field Bq > 10 -3 T. Thus, in such ultra-cold laboratory systems, a very weak external mag- 
netic field would make spin effects important for the formation of solitons, and the theory 
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FIG. 1: The Sagdeev potential ^{Z) (upper panel) and the profile of the solitary wave Z{£) (lower 
panel), for vq = 0.01 (dashed lines), vq = 0.5 (solid lines) and vo = 0.7 (dotted lines). The other 
parameters are e = 5, c s = 0.1, vb = 0.2 and \uj ce \ojc /^p e = 1- 




FIG. 2: The Sagdeev potential ^{Z) (upper panel) and the profile of the solitary wave Z{£) (lower 
panel), for e = 1 (dashed lines), e = 5 (solid lines) and e = 10 (dotted lines). The other parameters 
are vq = 0.5, c s = 0.1, vb = 0.2 and \u> ce \ujc /tOp e = 1. 
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FIG. 3: The time-dependent dynamics of the normalized magnetic field magnetic field b, for e = 5 
(left column) and e = 10 (right column). The other parameters are c s = 0.1, vb = 0.2 and 

\uJ ce \uJ C /^le = !• 



presented here could therefore be checked experimentally. 

In conclusion, we have investigated the effects of the quantum Bohm potential and the 
electron spin-1/2 on the existence of magnetosonic solitary waves in a magnetized quantum 
plasma. The solitary waves exist due to a balance between the nonlinearities and the dis- 
persion induced by the electron quantum diffraction/tunneling effects associated with the 
quantum Bohm potential. The spin introduces an additional negative pressure-like term in 
the quantum momentum equation, with the effect that solitary waves become wider and 
have shallower density depletions for larger values of the Zeeman energy e = [ibB /ksTe. 
We note that the spin term in the Sagdeev potential (13) can dominate the dynamics in the 
regime of C%, C\ <C C\v 2 B e. This regime corresponds to a dense quantum plasma with an 
ambient magnetic field, such that u) ce ujc "C oj pe and /ce(T e + Tj) <C ^bB$. Thus, the spin of 
the electrons collectively modifies the quantum dynamics of the MHD plasma significantly. 
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